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Abstract 

Influence of surface pair breaking, barrier transmission and phase difference 
on quasiparticle bound states in junctions with d-wave superconductors is 
examined. Based on the quasiclassical theory of superconductivity, an ap- 
proach is developed to handle interface bound states. It is shown in SIS' 
junctions that low energy bound states get their energies reduced by surface 
pair breaking, which can be taken into account by introducing an effective 
order parameter for each superconductor at the junction barrier. More inter- 
estingly, for the interface bound states near the continuous spectrum the effect 
of surface pair breaking may result in a splitting of the bound states. In the 
tunneling limit this can lead to a square root dependence of a nonequilibrium 
Josephson current on the barrier transmision, which means an enhancement 
as compared to the conventional critical current linear in the transmission. 
Reduced broadening of bound states in NIS junctions due to surface pair 

breaking is found. 

PACS numbers: 74.50. +r, 74.80.Fp 
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I. INTRODUCTION 



The important role in the Josephson effect of quasiparticle bound states locahzed 
in contacts between isotropic s-wave superconductors, is well known from studies of 
superconductor-normal metal-superconductor (SNS) junctions and quantum point contacts. 
Even for spatially constant order parameters of equal moduli on the two sides of a junc- 
tion (no surface pair breaking), Andreev reflection takes place on account of the phase shift 
between them. Current-carrying interface bound states with phase-dependent energies are 
formed at the junctionBHi'S. In short symmetric junctions, these bound states are known 
to carry all the Josephson current, while in asymmetric junctions the contribution from the 
continuous spectrum is of importance as well. 

The situation becomes more complicated in junctions involving d-wave superconductors 
where a quasiparticle, depending on its momentum, can see both substantial modulus dis- 
tortions or (and) a sign change in the order parameter in a reflection or a transmission 
process at a junction. Several kinds of bound states occur in this case, each with a different 
dependence on the quasiparticle momentum according to the incoming and outgoing quasi- 
paticle trajectories as well as the crystalline orientations on the two sides of the junction. 
In particular, low energy interface bound states are of interest, associated with changes of 
sign of the order parameter in reflection or transmission events. They become dispersionless 
zero-energy states both in the limiting case of zero transmission (impenetrable wall) or (and) 
in the opposite limit of a ballistic junctionii'iil. In certain conditions, zero-energy (or low- 
energy) bound states can result in an anomalous low temperature behavior of the Josephson 
critical currentiS and in characteristic peaks and jumps in the I-V curve JHI. In normal 
metal-insulator-superconductor (NIS) junctions they lead, at sufficiently low temperatures, 
to a zero-bias conductance peakiSiBS. 

As incoming quasiparticles with a momentum along the interface normal see the same 
d-wave order parameter as the outgoing ones, the bound states they occupy are to some 
extent analogous to those in junctions with isotropic s-wave superconductors. These bound 
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states can dominate the charge transport across a junction for certain orientations of d-wave 
superconductors if the transmission coefficient is sufficiently selective, limiting the transport 
of current to quasiparticles with momentum directions close to the interface normal, as is 
always the case for thick junctions. The difference between s-wave and d-wave cases, in 
general, is present even for the momentum orientation parallel to the interface normal being 
associated with surface pair breaking substantially more pronounced in anisotropically paired 
superconductors. As compared to isotropic s-wave superconductors, boundary conditions for 
anisotropic order parameters are very different because of surface pair breaking even within 
the Ginzburg-Landau theory001lilli. Depletion of the modulus of the order parameter in 
the vicinity of a contact modifies the bound state energies seen in non-selfconsistent models 
with constant absolute values of the order parameters. Besides, additional bound states can 
appear for quasiparticles in an effective potential well formed by the spatially dependent 
moduli of order parameters on the two sides of the junctionli§il. Since surface pair breaking 
depends upon crystal-to-interface orientations, it can strongly modify corresponding angular 
dependences of bound state energies and, in particular, the critical current obtained for a 
non-selfconsistent spatially constant order parameterJ^^. 

Assuming order parameters with spatially constant moduli on both sides of a junction, 
calculations of the energies of current-carrying quasiparticle bound states can be found in the 
literature admitting the presence of a phase shift and any value of the transmissionBUBIi. 
The effect of surface pair breaking combined with a phase shift and the influence of a flnite 
transmission on the interface bound states have, however, to the best of my knowledge, not 
been studied yet. Below, in Section II, I shall develop an analytical approach for studying 
the combined effects of surface pair breaking, a phase shift and the transmission coefficient 
on bound states, localized, in particular, at a contact between two (i-wave superconduc- 
tors. The approach is based on the quasiclassical theory of superconductivity focused on 
the problem of interface bound states. A key point of the consideration is that retarded 
propagators at the interface take quite large values (pole-like terms) at energies u close to 
an interface bound state energy esiPf)- Expanding propagators in powers of {u — SBipf)) 



one can introduce, in a first approximation, an ansatz for tlie bound states, wliicli reduces 
the Eilenberger equations and tlie normalization condition to one scalar differential equation 
for a quasiclassical phase rj to be completed with asymptotic and boundary conditions for its 



solutions. This equation was derived earlier in Ref. |T5| for the particular case of impenetrable 
boundary and real order parameter. Actually the equation has a more general character, 
as it is associated with the approach where the Eilenberger equations transform to a scalar 
Riccati equationiil. It has the same form as that obtained on the basis of the Bogoliubov-de 
Gennes equations within the WKB approximationilli, being applicable both to descrete and 
to continuous spectrum. Boundary conditions formulated below for the quasiclassical equa- 
tion are new. It is remarkable, that the boundary condition for the quasiclassical quantity 
7], in accordance as it is with Zaitsev's boundary conditions, can be formulated separately 
from all other quantities, just like the one-scalar boundary condition for the equation for 
7] mentioned above. It is a great deal simpler as compared to Zaitsev's original form and 
admits analytical results. The same boundary condition for rj can be derived also starting 
from the boundary conditions for the Andreev equations irrespective of whether a descrete 
or continuous spectrum is considered. 

On this basis I study interface bound states in both SIS' (see Section III) and NIS 
(see Section IV) junctions. It is shown for SIS' junctions that surface pair breaking results 
in reducing the energies of low energy bound states which can be taken into account by 
introducing an effective order parameter for each superconductor at the junction barrier. 
For interface bound states near the edge of the continuos spectrum, the effect of surface pair 
breaking turns out to be more interesting, resulting in splitting the bound state energies. 
In the tunneling limit this can lead to a square root dependence of the nonequilibrium 
Josephson current on the barrier transmision, which means an enhancement as compared to 
the conventional critical current linear in transmission. For NIS junctions the influence of 
surface pair breaking on broadening the bound state is considered. 
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II. QUASICLASSICAL THEORY OF INTERFACE BOUND STATES 



A. Ansatz for bound states 

Quasiclassical theory of superconductivity is based on Eilenberger's equations for the 
quasiclassical matrix propagator. In the case of a clean singlet anisotropically paired super- 
conductor the equations for the retarded propagator reduce to the following 2x2 matrix 
form: 



-Tin . 



(1) 
(2) 



Here, e, pf, and A are the quasiparticle energy, the momentum at the Fermi surface, 
the Fermi velocity and the order parameter matrix respectively. A "hat" indicates matrices 
in Nambu space and 7^3 is a Pauli matrix in this space. The propagator g and the order 
parameter matrix A have the form 
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/+ -9 



and 



A 



^ A^ 



V 



-A* 



(3) 



Henceforth the superscript R is dropped for simplicity. The boundary conditions for the 
quasiparticle propagators at a smooth interface with transmission D{pf) = 1 — R{pf) {R{pf) 
the reflectivity coefficient of the interface) are given by Zaitsev's relations (see Refs. |23|, ^4]) 
which can be written in the following matrix form 



Sl, Sr 



TT di 

2 



di — dr 



(4) 
(5) 



with a = (1 - R)/{1 + R), §1^,.) = 9i{r){PfMr)) + gmiPj^i^r))^ = 9i{Pf,i) - 9i{Pjj), dr = 
~ 9r{p f,i) + grip J i) , and the propagators are taken on the left or the right side of the interface. 
Equations (|^), (|) connect, at the interface, the propagators of an incoming quasiparticle 
from the left and the right sides of the interface with momenta p/ / , r and the propagators 



of the reflected quasiparticles with the momenta p^^ p^^. For specular reflection, the 
momentum parallel to the interface is conserved, i.e., = = ptj = ^1^,- For a 
complete determination of the quasiclassical propagator one has to take into account that 
deep inside the superconductor the propagator approaches its bulk value. 

In the presence of a quasiparticle bound state with the energy ^^(P/), manifesting dis- 
persion dependence on the Fermi momentum pj, the quasiclassical propagator g has a pole 
at e = esipf)- Following Ref. one can introduce the residue of the propagator g 

g{pf,r;eBiPf)) = 1™ [(£-eB(p/))^(p/,r;£)] , (6) 

which is finite, satisfies the same transport equation (P as g, but completed with the relation 

[l{pf,r;eB{Pf))f = , (7) 

rather than the normalization condition @. 

Linear boundary relations (|^) being applied to d, remain unchanged. At the same time 
nonlinear boundary conditions (^), taken at a bound state energy, simplify in a significant 
fashion. Terms containing multiplications of three propagators dominate in (^, if each 
propagator is well described by a large pole-like term. Other terms in Eq.(^, with only two 
propagators, can be neglected under conditions in question. Then Eq.(|D reduces to 

'■2 ft '~ 

5/ - ^{ShSr} = . (8) 

As the left hand side of matrix equation (j^) is proportional to the unit matrix, it leads to 
one independent scalar equation only. 

Eilenberger's equations for g can be solved in terms of the following ansatz: 

f+{pf,x;eB{Pf)) = g{pf,x;eB{Pf))exp{-ir]{pf,x)) , 
f{pf,x;eB{Pf)) = -g{Pf,x;£B{Pf))exp{iri{pf,x)) . (9) 

This ansatz was introduced earlier in Ref. ^ for the particular case of a real order parameter 
and impenetrable wall. In general, the substitution (||), satisfying Eq. (^, allows for the 



quantity 77 (as well as g{pf,x;eB{Pf)) and esiPf)) to take complex values r] = i]' + irj" 
{sBipf) = £'BiPf)~^'^^BiPf))- Complex values of r] and esipf) imply, in particular, broadened 
quasiparticle bound states (due to finite quasiparticle lifetime) discussed in the last section 
of the present article for NIS-j unctions on account of finite transmission. 

Introducing the phase of the complex order parameter A one 
obtains from ([l|), (|^) with substitution (||): 

g{pf,x;eB{Pf)) = ^o(P/, £b(p/)) exp J \^{Pf sin {r]{pf,x) - (j){p f , x)) dx 







(10) 

together with the following equation for t]: 

— ^^d^r]{pf,x) + 6 B{Pf) - \A{pf,x) \ cos {r]{pf,x) - (j){pf,x)) = . (11) 



According to Eq. (p!OD , under some conditions the residue g vanishes exponentially in the 
bulk of the superconductor. Then, the quantity g describes a quasiparticle state bound to 
the interface. Furthermore, equation (|n|) coincides with one obtained many years ago on 
the basis of WKB approximation for Bogoliubov-de Gennes equationsili, being applicable 
both to descrete and to continuous spectrum. For continuous spectrum rj{pf,x) — (j){pf,x) 
should be a purely imaginary quantity in accordance with Eq. ([lO|) . The same equation (pUf) 
is known to appear as well within the approach transforming the Eilenberger equations to 
a scalar Riccati equation^. 

The asymptotic condition for rj, which garantees solution (0) to vanish in the bulk, 
takes in the right half space {x — > +00) the form 

v},x(P/,r) sin (<oo(P/,r-) " 0r,oo(P/,r)) > , (12) 

while in the limit x ^ —oo 

V/,x(P/,«) sin {v'i,ooiPf,i) - 0z,oo(P/,O) < . (13) 
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Substitution @ essentially simplifies relations 
following boundary condition for equation (|11|) 



resulting, in particular, in the 



Rsin 



'vi,oiPf,i) -VifiiP 



sm 



'Vr,o{Pf,r) - VrfiiPj^^.. 



' Vl,0{Pf,l) - Vr,0{Pf 



D sin I — '■ '■ sin ' — — 



VlfiiPf,) -^r,o(P/,, 



(14) 

This boundary condition holds for any value of the transmission coefficient. In the tunneling 
limit, ^ 1, it is convenient to transform Eq.(|l4D to the following equivalent relation 

'%o(Pf,l) -Vrfi{Pf,y 



n • ( Vl,oiPf,l) - Vr,oiPf,r)\ . \ ■l'^''''flf,l> ■rr,v\t^f J 

U sm — sm ' 



sm 



'viAPf,i) - viAPf,i) 



'Vr,o{Pf,r) - VrfiiPj^,., 



sm 



(15) 

It is remarkable that the boundary condition for the quantity rj can be formulated separately 
from other quantities, simply as a boundary condition for equation (pHI). Other boundary 
relations for quantities entering ansatz (P) are given in Appendix. 



Relation ([1^) can also be derived within a more general framework independent of 
whether a descrete or continuous spectrum is considered. For this purpose one can rep- 
resent the Andreev amplitudes in the form 



u{pf,x) 



\ ( gi»7(p/,2:)/2 \ 



V 



J(,{pf,x) 



(16) 



where ri{pf,x) and ^{pf,x) are, in general, complex. Substituting ( p^D into the boundary 
conditions for Andreev amplitudesii, one obtains Eq. (|T^) as a separate boundary condition 
for ri{pf,x)^. In the particular case of descrete spectrum ri(pf,x) is a real quantity, while 
C,{pf,x) is purely imaginary leading to exponentially decaying asymptotic behavior of the 
Andreev amplitudes. Furthermore, functions entering the expression for the quasiclassical 



matrix Green function and satisfying the Riccati equation can, in general, be represented 

Boundary conditions for these functions are evidently 
directly related with Eq. ([T^). It is worth noting, in addition, that Eq.(pJ]) transforms to 
the Riccati form by introducing the new function (3 = ta.n{ri/2). 

If ri{pf,x) is a solution of Eqs. (|TT|), ( P^ with the energy e{pf), then tt — ri{pf,x) is a 
solution of the same equations with ~e{pf) for the system with a given pf and complex 
conjugated order parameter. So, for a given pf quasiparticle descrete spectrum of a system 
can be, generally speaking, asymmetric with respect to the Fermi surface under the condition 
0^0. 



B. Positions of poles and residue values 

As it is shown in this section, energies of bound states, entering denominators of pole-like 
terms in the expressions for the propagators, on the one hand, and residues of those pole-like 
terms, on the other hand, can be expressed via ri[pf,x) for a given A{pf,x). 

Since dxi]{pf,x) vanishes in the bulk {x —>■ ±00), one immediately gets from Eqs. (pUf ) 
the relation between the bound state energy esipf) and the quantities rjodPf), Aoo(p/) in 
the bulk of the superconductor: 

eB{Pf,r) = |A^(p/,^)| COS {rir,oo{Pf,r) " 0r,oo(P/,r)) = \Al^{p^^^)\ COS (r?r,oo(p^_J " <Pr,oo{Pj^^)) 
= COs(r/;,oo(P/,/) - 0/,oo(P/,O) = \^LiPf^i)\ COS {vi,ooiPf,i) ' ^t^lAPf)) ■ 

(17) 

As a consequence, bound states might exist for a given momentum direction only below the 
band edges for the momenta P/,z(r) and Pjn^^y i-e-; for 

\eB{Pf)\ < min {|A;^(p^,)|, |A^^(p^,,)l, \^L{Pj,)l • (18) 

Furthermore, for a frequency near the bound state energy esiPf), the quasiclassical 
Green's functions can be expanded in powers of {u — esiPf))- Taking into account ansatz 
(I), one has 
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9iPf,x;uj) = ^^^^l^^^j^^^ , (19) 



[uj-eB{Pf)+iS) „=o 



n=0 



riPf. x; ^) = + ^ ^^^"^(P/' ^' ^b(p/))(c. - .i.(P/))" . (21) 

Further I introduce the quantities 

h{pf,x-u) = 1 {r{pf,x;uj)e'^^Pf'^^ ± /(p^, x; u;)e-'^(P/'^)) . (22) 

According to (|T9D -(pTD, f^{pf,x;uj) has no singular part (no pole-hke term), while the 
singular part of f-{Pf, x; u) coincides with the one for g{pf, x; u), that is 

f~{Pf,x;eB{Pf)) = g{Pf,x;eB{Pf)) , f+{pf,x;eB{Pf)) = 0. (23) 

In addition, substituting expansions (0)-(^) into the normalization condition (0) and 
equating terms inversly proportional to (tu — £b(P/)), one finds 

f^^\pj,x-eB{Pf)) = g^'\pf,x-eB{Pf)) • (24) 

Taking into account (]22|)-(p4D one derives at the following relationship (for x > 0) from 
the Eilenberger equations 

CO 

ff\pf,x;eB{Pf)) = — dxg{pf,x;eB{Pf)) + f+^oo{Pf;£B{Pf)) ■ (25) 

Vx J 

X 

The bulk value of function ff\pf, x; SBiPf)) can be easily found: 

/SL(P/; SBiPf)) = __JI}^^M== sin iv^ip,) - 0^(p,)) . (26) 

|Aoo(P/)P -4(P/) 



Substituting Eqs.(|ToD, (|26|) into (||) I get 
ff{Pf,x = 0;£b(p/)) = i sgn{vx) \A-\pf,0)\ goiPf^SBiPf))- 



^^|Aoo(P/)| sin(r/^(p^)-0«.(p;)) , (27) 



Aoo(p/)P-4(P/) 
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where 



1 

|A(P/,0)| " 



2 



y"exp(-- — - — -y"|A(p/,X2)|sin(?7(p/,a;2) -0(p/,a;2))rfx2 Ma;i . (28) 



\^fAVf)\{ \ ^fAVs) ^ 



In terms of Eg . (|27|) , (Pq), one can easily obtain expression for the residue of the propa- 
gator taken at the impenetrable wall. According to the boundary conditions jor an impen- 
etrable wall, quantities f^\pf,x = 0;eB{Pf)) and g{pf,0;eB{Pf)) = %{p f ^ ^ b{p f)) ■, taken 
for incoming momenta are equal to the same quantities of the outgoing ones. Then one 
obtains from (0), (p^ ) and (|T2|) the following expression for the residue ^o(P/) of 
the propagator g{pf,x > 0]u) taken at the wall (x=0) and at the bound state energy: 

27r|A(p„0)||A(p 0)1 

9o{Pf,£BiPf)) = ~ . (29) 

|A(p^,0)| + |A(p/,0)| 

The positive sign of the residue stipulates a positive contribution from each bound state 
to the angle resolved local density of states at the wall. In accordance with Eqs. (P^, (P^D, 
the residue is fully determined by the quantities A{pf,x), ri{pf,x). The quantity rj{pf,x) 
obeys differential equation (|ll|), while the position dependent order parameter A(p/,a;) 
needs to be determined self-consistently after specifying a particular form of the pairing 
potential. The problem of self-consistent spatial dependence of (i-wave order parameter near 
the interface was explicitly studied numerically, for example, inlllilii. The self-consistent 
space dependent order parameter A(pj, x) is considered to be given throughout this article. 

For the particular case of midgap states and real order parameter one haslll esipf) = 0, 
sgn[Aoo(p/)] = -sgn[Aoo(p^)], v{Pf,x) = r]oo{Pf) = f sgn [v/,^(p/)] + 0oo(P/), 0(P/,a;) = 
0oo(P/) = ■J (1 ~ sgn[Aoo(p/)]). Expressions (p8|), (^) then reduce to those obtained earlier 
in Ref. |TT|. 
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III. BOUND STATES IN SIS' JUNCTIONS 



In order to calculate a bound state energy within the framework of the approach devel- 
oped in the preceding section, one should first find solutions to equation (pA]), which satisfy 
both asymptotic relations (0), (|13|) and the boundary condition (|1^). The bound state 
energy esipf) is directly associated with the asymptotic value rjaoipf) according to Eq.(|T^). 

Explicit analytical expressions for the bound state energies esiPf) will be found below 
under certain conditions in terms of spatially dependent order parameters on both sides of the 
interface. As was demonstrated in Ref. for the particular case of a massive supeconductor 
confined by impenetrable wall, effects of surface pair breaking for various (although not for 
all) quasiparticle trajectories result only in weak deviations of the space dependent quantity 
r]{pf,x) from its asymptotic value rjaoipf)'- \5'iliPfi^)\ = IviPfi^) ~^oo(P/)| ^ 1- In the 
case of finite transparency of a junction of two halfspaces one can assume this condition as 
well, linearize equation ([TT|) with respect to 6ri{pf,x) and find the solutions of Eq.(|ll]) on 
account of asymptotic conditions ([1^), ([T3|) and Eq.(|l^: 

Vi{Pf,i,x) = MPli) - sgn(v^,^(p/,/)) arccos ( ) + / dx^x 

/ \A^{Pf^i,xi 



A\p„,X2)\dx2\ , (30) 



Vr{Pf,r, x) = (pr{Pf,r) + sgn(v^^^, (p/,^) ) arccos 



£B{Pf,l) \ (p/,. 



x(P/,r) 



7: 



n+oo 

/ dxix 

J X 



\A''{pf^r,X2)\dX2 



(31) 



UPf,r)\\ \^loiPf,r)\ 

Effects of supercurrent flowing across the junction (along the x axis) can be taken 
into account by adding the spatially depending term 2mVsX into the phases (f>i{r){Pf,i{r)) in 
(pTp. Then the corresponding solutions are obtained from (0), (|3TD after the substitution 

(pl{r){Pf,l{r)) + '2mVsX. 



£B{Pf,l) SB{Pf,l) - vf'\pf,lir))mVs, 0«(r)(P/,/{r)) 
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Bound state energies can be found now combining solutions (|30D, ( pi]) and the boundary 
condition (pISl). The equation obtained within the approach should be linearized with respect 
to the effects of surface pair breaking ( discribed by spatial integrals in (pO|), (pl|)) both in 
the right and the left superconductors. 



A. Low energy bound states 

Let us consider bound states with low energies 

esiPu) « min{|AL(p/,0|, |AL(p^^,)|, |A^^(p;,,)|, |A^^(p^^J|} . 

Linearizing Eqs.(pO|), ( pT]) with respect to the small parameters eB{Pf,i) /\^''^\pf,i{r))\ 
and taking into account the presence of a supercurrent, one gets after simple transformations: 

ViiPLi, = MPf,i) - lTSgn(v^,,(p/,z)) + ^ {- , (32) 

2 \A{pf^i,x)\ 



I ^ [^b{Pu) - y^fiPf,r)mVs) sgn(v^_,(p/,,)) 

r]r{Pf,r, x) = MPf,r) + 7TSgn(vy_^(p/,,)) - ^ , (33) 

2 |A'^(p/,^,x)| 



where the following quantities are introduced 

1 2 ( 2 



1 2 f 2 \ 

T77 ^ = T^f — \\ / ~T^f — \\ / \^'^{Pf,X2)\dx2 dxi . (35) 



I first assume that m the limit D ^ midgap states arise on both sides of the harrier 
plane for given crystal to surface orientations and quasiparticle trajectories considered. For 
simplicity, let the phases of the order parameters (j)i{pf^i), (p^ipf^r) in both superconducors 
be spatially constant in the absence of a supercurrent. 

Then, allowing for the difference between the phases of right and left comlex order 
parameters, these phases may be written under the conditions considered as follows 

(f)l{pf^l) = (j)l + 2mVsX , (f)r{Pf,r) = 4>r + T^Sii +2mVsX , 

(36) 
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where two different cases i = 1,2 are taken into account. The superffuid velocity Vs is 
assumed to be positive (negative) for the supercurrent flowing along (opposite to) the x- 
axis. 

Expanding boundary condition ([T5| ) with the substitution (^2]), ( p3D with respect to 



{eB{pf,i) — T^f^s^''f^\pf,i{r))^, results in a simple quadratic equation for the bound state 
energies: 

, JA'(p/,,0)| + |A'(p^^,,0)| , , , ( \yUPf,i)\ 1^/-^^^/.^^ 

^B{Pf,r)—jr^ nMiXU ^;^sgn(v^,^(p/,i)) - mt;, ' 



|A'(p/,.,0)||A'(p^^^,0)| w,--/.- ^^|A'(p^,,0)| |A'(P^,,0)| 

, jA-(p/,.,o)| + |A-(p,,,o)| ( \viM,r)\ I-Upj;. 

|A'^-(P..0)||A^(P,^,0)| (^.(^ - 

= 4Dsing + ^5,)cos(^-^5.) . (37) 

For the sake of simplicity let D ^ 1 in Eg . (|37|) . 

Two particular solutions of Eg . (1371) are of special interest. For a "symmetric" tunnel 
junction (STJ) identical superconductors with the same orientations are situated on its 
left and right sides. Subscript i = 1 in Eqs.(^), (^) corresponds to this case. Absolute 
values of order parameter on both sides of the STJ taken at the same distance from the 
interface, are equal to each other, even in the presence of a phase difference and supercur- 
rent: |A'(pf;,— x)| = |A''(p ,x)|, |A'(p ,— x)| = I A''(pf x) |. Then analogous equali- 
ties follow from Eqs.(|3^), (p5|) for effective order parameters: \K\pf^i, —x)\ = |A^'(p^^, a;)|, 
l^KPfP ~^)\ = |A'(P/,r5 a^)|- Under these conditions one gets 



2|A-(p,,.0)||A-(p,^^.0)| 
" |A'(p,.„0)|+|A'-(p,..,0)| 



COS — 

2 



+ 



v;^,,(p;,.)|A-(p^^^,0)| + v-,,(p^^J|A-(p;.„0)| 



mvs ■ (38) 



|A^(p/,.,0)| + |A''(p^^,,0)| 

For a "mirror" tunnel junction (MTJ) i = 2 orientations of identical superconductors 
can be obtained from each other, by definition, by a reflection with respect to the junction 
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barrier plane. Then \A\pf^i, —x)\ = \A'^{pf^r,x)\, \A\p , —x)\ = \A'^{p ,x)\ and one 
obtains the corresponding values of bound state energy from Eqs.(p4D, (p5D, (pT]): 



W{Pf^rMW{Pfy^)\ 



|A^(p;,.,0)| + |A^(p^_^,0)|\ 



W sin^ 



+ 



(39) 



As Vg can depend upon e^, relations (|38|) , (p9|) still represent, generally speaking, imphcit 
equations for eb- However, factors in front of mvs in Eqs.(p8D, ( pQ]) can be extremely small 



or even vanish (as in the case | A''(j?^^, 0)| = |A''(pj^r, 0)|, v'f^^i.Pf^ 



(Pf^r))- Then 



one can disregard the influence of the superfiow on the bound state energy, and get fairly 
simple results from Eqs.(p3), (pQ]) 



^B[Pf,i 



2\A^{pj,r,0)\\A^ip 



\A^{Pf,,M + W{Pf, 



D I 



cos 



sm 



for ST J, 



for MTJ . 



(40) 



One should note the relation Eb oc VD and the essentially different dependences of the 
bound state energies upon the phase difference (j) for ST J and MTJ. According to (|38D-(P0[), 
zero energy bound states at an impenetrable wall {D 0) in the absence of supercurrent, 
split into two low energy levels on account of the effects of nonzero (low) transmission and 
phase difference. For STJ the bound state energy has its maximum for = and the 
principal reason for the bound state energy shifting away from the midgap position is the 
finite transmission. For MTJ the bound state energy differs from zero for nonzero (p taking 
its maximum value at = vr. This difference is evidently related to the additional phase 
shift TT acquired by paired quasiparticles crossing the junction and seeing different signs of 
the gap functions for the given momentum direction in two superconductors in MTJ. Let 
now midgap states take place in the limit D — > only on one side of the barrier plane, 
for instance, in the right superconductor for quasiparticle trajectories considered. Then 



15 



expressions for the phases of complex order parameters differ from Eq.(^6]) and have the 
following form 



MPf,i) 



+ 2mVsX 
+ 2mVsX 



iPf,r] 



+ nSii + 2mVsX , 
— 7r5j2 + 2mVsX . 



(41) 



In the case of an impenetrable wall there is no solution to equation (|TTp in the left 
half space x < for sufficiently small energy. Such a solution arises, however, at finite 
transmission, being induced by the proximity effect on account of the corresponding solution 
for the right half space. These solutions have the same form (pOf ), (pl|) with new relations 



Boundary condition (|T^) then reduces to 



7r 



Vrfi[Pf,^ 



(P/,rO - ^^M'^lxiPf'r)) 



TC 



Vr,0[P 



■f,r> 



-ip^Dsincf) . (42) 



In terms of Eqs. (^2]) , (|53D-(p8D, one can justify that, indeed, gi^ oc -D^r,o both for the 
incoming and outgoing momentum directions. 

Substituting r]rfi{Pf,r), Vr,o{Pr ) from Eq.(|^) into Eq.(^, one gets the following expres- 
sion for the bound state energy 

|A^(p/,.,0)||A^(p^^^,0)|sgn(v^,,(p^^;; 



|A^'-(P/,.,0)| + |A^Xp^,,,0)| 



-1 



D sin I 



UPf, 



(43) 



|A''(P/,.,0)| |A^'(P^,,,0)|^ 
Results obtained in this section show that the effect of surface pair breaking on the low- 
energy bound states, can be taken into account by introducing the effective surface values 
of order parameters definded in Eq.(^), (^). If one disregards for the time being the 
surface pair breaking at all, considering moduli of the order parameters being independent 
of spatial coordinates near interfaces, then it follows from Eqs. (|3^) , (^) |A(p/,x)| = 



|^oo(p/)| = const (x). For crystalline orientations, when |Aoo(p/,r 
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l^oo(p^ ,,)|, one gets 



from Eq.(^ in this particular case e^'^{pf^r) = ±|Aoo(p/,r)| 



cos ■ 



c-MTJ 



sm 



. According to (0), p5|), the relation \A{pf,x)\ < |Aoo(p/)| holds 
as a consequence of surface pair breaking. Disregarding this effect leads to an overestimation 
of the shift of the zero energy bound states brought about by the finite transmission of the 
junction barrier. Such an overestimation can be quite noticeable. Moreover, for smooth 
surfaces the suppression of the order parameters depends on crystal to surface orientation. 



The effective surface quantities (|3J), ( p5|) can manifest qualitatively different dependences 
upon crystal to interface orientation as compared to the dependences of the order parameters 
in the bulk. 



B. Splitting from surface pair breaking of the bound states near the edge of the 

continuous spectrum 

In junctions of isotropic s-wave superconductors, the order parameters with incoming and 
reflected quasiparticle momenta are identical. In the s-wave case one can usually disregard 
surface pair breaking and consider the order parameter as spatially constant up to the 
interface. Then the energies of the interface bound states in symmetric junctions, as it is 
well known, are as foUowHi: 



esie) = ±|A| ^1 - D{e) sin^(0/2) . (44) 
It is explicitly indicated here that transmission depends upon the angle 6 between momentum 



direction and the interface normal. For tunnel junctions {D <^ 1) energy (|4^ ) lies close to 
the edge of continuous spectrum. 

In preceding section the presence of quasiparticle trajectories was assumed, where in- 
coming and outgoing quasiparticles see the order parameter values of opposite signs. The 
corresponding states have quite low energies in the tunneling limit. Interface bound states 
in the vicinity of continuous spectrum, however, always arise in tunnel junctions between d- 
wave superconductors as well. This is the case for quasiparticle trajectories with momentum 
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directions along the interface normal, where incoming and reflected quasiparticles always see 
the same order parameter of an anisotropic singlet superconductor. Interface bound states in 
d-wave superconductors for trajectories in close vicinity of the interface normal are to some 
extent analogous to the ones in isotropic s-wave case. For special crystalline orientations 
such as 45° in dr^2_y2 superconductors all those trajectories reduce to the only one along the 
interface normal. 

If one considers a contribution of bound states to the current across the junction, then 
the particular dependence of a transmission upon 9 turns out to be especially important. 
For the transmission gradually diminishing with increasing 6', like in case of sufficiently thin 
barriers, trajectories in a wide range of 9 can substantially contribute to the current. By 
contrast, if the transmission noticeably differs from zero only for small 9 (this is the case for 
thick barriers), the contribution from trajectories with sufficiently small 9 entirely governs 
the current. Below bound states injunctions between cZ-wave superconductors are considered 
just for the trajectory along the interface normal. 

Bound states are formed due to Andreev reflection processes resulting from a phase 



difference in the order parameters of a symmetric junction. According to (Q), the bound 
states merge into the continuum in the limit D — 0. However, this is not actually the case as 
a result of a suppression of the moduli of the order parameters in the vicinity of the barrier 
plane. For a smooth impenetrable wall, the bound states still exist on account of surface 
pair breaking for momentum direction along (or sufficiently close to) the surface normal for 
almost all crystal to surface orientations of a rf-wave superconductorEl. Bound states of this 
kind arise even for a weak suppression of the order parameter. Then, however, there must 
be a tiny distance from the energy of a bound state to the continuum. For isotropic s-wave 
superconductors this applies to any quasiparticle trajectory. 

Inhomogeneous modulus of the order parameter as such near an impenetrable wall can 
be considered an effective potential well for quasiparticles. Finite transmission across the 
interface of a junction, as opposed to the impenetrable wall, results in a double well structure 
and allows for the interplay of the phase difference and the inhomogeneous moduli of the 



order parameters in forming interface bound states. It is shown below, that sphtting of 
bound state levels takes place in the double well structure for the quasiparticle trajectory 
along the surface normal, at least in the tunneling limit. 

One can easily see from (p!?]) for positive bound state energies near the continuum, that 
in "symmetric" or "mirror" tunnel junctions all asymptotic values {rjodPf) — 0oo(P/)) are 
small. Assuming in addition \ri{pf^x) — ri^{pf)\ <^ 1, one can get from Eq.(|ll|) in the first 
approximation the following expressions for the corresponding interface values ri{pf,x = 
0) = VoiPf) ■■ 



Vi,o{Pf,i) = -sgn(v^_^.(p/,/)) 



+ 



|A-oo(P/,; 



dx 



(45) 



VrfiiPf,i) = + sgn(v}^^(p/,/)). 



\^'{Pf,r,x) 
|ASo(P/,r)l 



dx . 
(46) 



Substituting these expressions into the boundary condition (p!5| ) and linearizing in 
(voiPf) — Vo{Pj)^i one obtains for STJ and MTJ the following bound state energies: 



^B[Pf,r) 

ial(p/,.: 



1 



I A r / M I dx ±VDsm\ - 

\^ooiPf,r)\ J V2, 



(47) 



The minus sign in front of the VD - term in (^7]) is admissible only when the modulus 
of this term is less than that of the first term in the square brackets. This implies a 
relatively large role of interface pair breaking in forming Andreev bound states as opposed 
to the contribution from the phase difference. As can be seen in Eq. (|17|) , under these 
circumstances the combined effect of the phase difference and the interface pair breaking 
results in a splitting of the bound state energy. In the opposite limit of absence of noticeable 
surface pair breaking, when \ A{pf, x) \ = |Aoo(p/)|, the first term within the square brackets 
in Eq. (|47|) vanishes. Then Eq . (|47|) reduces to the conventional result Eq.(|^). On the 
other hand, for an impenetrable barrier Eq . (^71) takes the form found in Ref. If order 
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parameters have equal moduli and opposite signs on the two sides of the junction, then the 
bound state energies are described by Eg. (|i7|) after replacing sin by cos (^^^■ 

The splitting considered above is formally analogous to the splitting discussed in Ref. ^ 
for long SIS junctions and specially constructed SNS junctions with double- well or double- 
barrier structures. However, the physical reasons for the splitting and the conditions for 
its observations considered in this section, drastically differ from what was studied earlier, 
where the effect of interface pair breaking was always neglected. The double well structure 
discussed in this section is inherent in junctions with noticeable interface pair breaking. One 
should note that the first term in square brackets in Eq.(^) is considerably less than unity 
both for s-wave superconductors as well as for d-wave superconductors. In the former case 
this is entirely due to a small suppression of the order parameter at the interface. The d-wave 
order parameter taken for a momentum along the surface normal, is small in itself for those 
crystal to interface orientations, which entail substantial surface pair breaking. However, at 
least for the d-wave case the first term in square brackets in Eq.(^7Dcan exceed the second 
one, allowing an observable splitting (the estimations follow from Fig. 3 in^l). 

The supercurrent fiowing via individual bound state (|T7| ) can be written as 

UP K Jo |AL(P..) 



I± oc e|A^(pj,,)|V/^cos ( - 



dx 



(48) 



For sufficiently small transmission the second term in the square brackets dominates, re- 
sulting in a supercurrent proportional to As it is known for resonance Josephson 
coupling in double well structures@, in equilibrium such anomalous terms coming collec- 
tively from all bound and extended states cancel in the expression for the total current. 
For nonequilibrium quasiparticle occupation of bound and (or) extended states, however, 
one can expect observable manifestations of the anomaly discussed above. As compared to 
the conventional behavior linear in the transmission coefficient, a V^D-term results in a in- 
hancement of nonequilibrium Josephson current from surface pair breaking in the tunneling 
limit. For d-wave superconductors, the broadening of bound states on interface roughness 
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and impuritiesEJ is of importance, since can smearing out the splitting. 

It is not entirely clear whether the A/D-behavior of the critical current vanishes in equi- 
librium within more general assumptions than made above. Such a study is in progress 
now. 



IV. BROADEN BOUND STATES IN NIS JUNCTIONS 

Andreev bound states, localized near an impenetrable wall in a d-wave superconductor, 
become slightly broadened in the case of small but finite barrier transmission of an NIS 
tunnel junction. In this case, the retarded propagators have no singularities at any real 
value of the energy , since positions of poles move from the real axis in the complex energy 
plane. This shift is linear in the barrier transmission coefficient, which is assumed to be 
sufficiently small. A relatively small imaginary part of a pole position is associated with 
finite quasiparticle life time for the quasistationary bound state. 

The other important feature of NIS (normal metal - insulator - superconductor) junctions 
is that the solution for the normal metal halfspace cannot be described by ansantz At 
the same time the dominating terms in propagators still satisfy Eq. (H) in the superconducting 
half space in the case of sufficiently small transparency of the barrier. 

As it follows from Eilenberger's equations, the diagonal component of the quasiclassical 
Green's function is constant throughout the normal metal up to the interface {g = —iir). At 
the same time superconducting correlations are present there due to the proximity effect and 
known to decrease exponentially toward the bulk normal metal (x < 0). Since imaginary 
part of a pole position should be negative for the retarded propagator, one can write for the 
anomalous Green's functions for x < 0: 

2eBiPf,i)x\ 



fNiPf,hX,eBiPf,i)) = fmiPf,i,SBiPf,i))ew [i- i , . i ■ 
Jn {Pf^pX,^BiPf,i)) = fm {Pf^p^BiPf,i)) exp j _^?£^ipii)^ 



(49) 



VfAP 



while {pj:pX,eB{Pf,i)) = 0, {pf^i,x,eB{Pf,i)) = . Here and below normal to the 
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interface velocity component v^j^^{pf^i) ( v^^^^ij)^^ ) is chosen to be positive ( negative ) for 
an incoming ( outgoing ) quasiparticle in the left halfspace. 

Evidently, propagators (^9|) in the left halfspace do not satisfy relation (^. There is, 
however, another important relation, which strictly holds in the normal metal halfspace 
and substantially simplificates Zaitsev's boundary conditions at the NIS interface: diagonal 
components of matrix di are equal to zero in the case considered. According to (H), diagonal 
components of vanish in this case as well. Thus, applying ansatz to the right halfspace 
and taking into account (|^), (^), one gets 



fm{Pf,i,eB{Pf,i)) = gr,o{Pf,r,eB{Pf,i)) 



JVrfiiPf,r) _ ^^Vr,o{Pf^r' 



fmiPf,r^B{Pf,l)) = 9r,o{Pf,r,£B{Pf,l)) 



(50) 



With relations just mentioned above the matrix boundary condition reduces only to a 
single scalar equation. In the particular case of this section one should not forget either that 
equations (^, (^) are valid only for sufficiently small a. For these small values of a the es- 



timation rfrfiipf^r) — VrfiiPjj.) oc « holds and the quantity agrfi{pf,r, £BiPf,i)) is sufficiently 
small and vanishes in the limit a — > 0. This can be seen, in particular, from (0), since in the 
limit of impenetrable boundary amplitudes fm{.Pf,i,£B{Pf,i)) , fNoi.Pfv^B{pf,i)) should van- 
ish. For sufficiently small transmission coefficient and therefore small ?7r,o(P/,r) ~ Vr,o{Pj^) 
a remarkably simple boundary condition follows from (H) 



Vr,o{Pf,r) - Vr,o{PfJ = -2za ^ -iD{pf) 



(51) 



The complex pole eB{pf,i) can be found now on the basis of equations (|ll]), (p^), (^TJ). 
For broadened zero energy bound states, one should assume 0j.(P/,r) = 0, 4>r{Pj^) = tt 
ensuring the existence of that state in the limit D ^ for the trajectory. Then the solution 
of Eq.(p!TD for the superconductor coincides with ( |33D in the absence of a supercurrent 
across the junction: 6rir{pf^r,x) = —eB{Pf,r)sgn(Y'j^^{pf^r)^/\A{pf^r,x)\. The definition 
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of \A{pf^r, x)\ = \A'^{pf^r,x)\ is giveii in (^). Substituting this solution into boundary 
condition (|5TD one gets a purely imaginary value eB{Pf,r) = ~i^{Pf,r) for the pole : 

vf ^ nf , |A(p,,,0)||A(p 0)1 

According to ([5^), the effects of surface pair breaking on the broadening of the zero energy 
bound state is taken into account by introducing effective order parameter values (^) taken 
at the interface. This effective order parameter reduces into order parameter of the super- 
conductor when neglecting surface pair breaking. This is entirely analogous to what was 
obtained in Sec. III. A for the effects of surface pair breaking on the energies of low energy 
bound states. In the particular case of no surface pair breaking Eg. (^2]) coincides with the 
result obtained in Ref. ^for crystal orientations when \A{pf^r)\ = |A(Pj^)|- 
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APPENDIX 

One can derive from the boundary conditions (||) with the substitution (j^, the 
following relations among interface values of the propagator g taken both for incident and 
(or) for outgoing momenta: 

gi,o{Pf,i) sm ( ) sm I 

'VlfiiPfl) - Vr,o{Pf,r)\ I VlfiiPfl) - VrfiiPfr^ 

9i,o[Pf^i) sm I — sm ' 



(53) 



g.,ip,.) sin hrAPf.)-mM] (^rAPf.) - V,oiP,,) 



'Vr,oiPfr) - Vl,o{Pf,l)\ ( Vr,o{Pfr) " ^'.o(P„y 

9r,0[Pf,r) I n 



(54) 



~ r ^ • (vi,o{Pf,i)-Vr,o{Pf,r)\ . (vi,o{Pf,i)-Vi,o{Pj^ij 
9i,o{Pf,i) sm sm ' 



'Vr,o{Pfr) - ^l,o{Pfi)\ I Vr,o{Pfr) ' Vr,o{Pf,^ 

9r,o[Pj^^) sm I — I sm ' 



(55) 



^.0(P/.)sin hr,iP,r)-V,.iPf,)\ (Vr,iP,r)-Vr,ip,.. 
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9iAPf,i) sin 



sm 



(56) 



^ ^ ^ . 'Vl,o{Pf,l) -Vr,o{Pfr)\ . (Vl,o{Pf,l) -Vl,o{PfiJ 

gifl{Pf,i) sm I ^ I sm ' ^ 



. , , . ,'Vr,oiPf,r) -VlfiiPfij 
-9r,o{Pf,r) I — I sm 



'Vr,Q{Pf,r) - Vr,0{Pf,, 



(57) 



'VlfiiPfi) - Vr,o{Pf,r)\ . I VlfiiPfi) - r]lfi{Pf,l)' 

9i,o{Pf^i) sm I — z I sm ' 



-^r,o(p_^^^)sin 



Vr,o{Pf,r)-'^lAPf,l)\ . ( VrAPf,r)-'^r,o{Pf,ry 



sm 



(58) 
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